Anderson localization of partially-incoherent light 
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We study Anderson localization and propagation of partially-spatially incoherent wavepackets 
in linear disordered potentials, motivated by the insight that interference phenomena resulting 
from multiple scattering are affected by the coherence of the waves. We find that localization is 
delayed by incoherence: the more incoherent the waves are, the longer they diffusively spread while 
propagating in the medium. However, if all the eigenmodes of the system are exponentially localized 
(as in one- and two-dimensional disordered systems), any partially-incoherent wavepacket eventually 
exhibits localization with exponentially-decaying tails, after sufficiently long propagation distances. 
Interestingly, we find that the asymptotic behavior of the incoherent beam is similar to that of a 
single instantaneous coherent realization of the beam. 

PACS numbers: 42.25.Dd,42.25.Kb,72.15.Rn 
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The phenomenon of Anderson localization was con- 
ceived in the context of disordered electronic systems 
however, localization phenomena have been exten- 
sively studied also in other systems [2l-|l7j including op- 
tics 0, H, [sl-flij and ultracold quantum gases p^Hl7| . 
For a direct observation of the phenomena, optical and 
ultracold atomic systems have some profound advan- 
tages over the condensed matter systems: the influence 
of the environment such as thermal fluctuations and 
phonons can be minimized to become negligible, and 
nonlinearity /interactions can be controlled and virtually 
turned off. Nonlinearity in optics can be controlled by 
the intensity of light, whereas interactions of ultracold 
gases can be tuned via Feshbach resonances. In con- 
trast, electron-electron and electron-phonon interactions 
are always present in condensed matter systems. These 
influences arc important since they affect interference, 
and Anderson localization arises from interference among 
multiple scattering events from disorder in the medium. 

Here we address Anderson localization of waves with 
imperfect coherence. We investigate whether an initial 
finite-size partially-incoherent wavepacket would spread 
through a linear disordered potential, or would its spread- 
ing be stopped by virtue of disorder? We demonstrate 
our findings on an optical (1+1)D potential, and discuss 
the implications of our results to other wave systems in 
nature. A superficial answer to the question above may 
be that, since incoherence destroys interference effects, 
a sufficiently incoherent wavepacket will diffuse through 
the random medium without ever being localized. How- 
ever, from the theory of coherence [Tsl it is known that an 
incoherent wave can be thought of as a superposition of 
coherent modes with stochastically varying coefficients. 
Each of these coherent modes is expected to undergo lo- 
calization, hence the entire wavepacket should localize. 
Our study reveals that coherence indeed affects the prop- 
erties of evolving wavepackets, in a sense that more inco- 
herent beams spread more while propagating through the 
random medium. However, if the system's eigenmodes 
are all localized, as is the case for all one-dimensional 



[(1-f 1)D] and two-dimensional [(2-|-l)D] fully disordered 
systems [l^, the incoherent wavepacket will eventually 
become Anderson localized. Finally, we find that a typi- 
cal instantaneous coherent speckled realization of the in- 
coherent beam exhibits the same asymptotic behavior as 
the time- averaged incoherent wavepacket. 

The idea that localization could be observed in optics 
dates back to the beginning of the 80s 0, Q . The experi- 
ments on so-called weak localization Q , which can be pic- 
tured in terms of a coherent backscattering process, were 
soon to follow Q . Experiments on strong localization 
in random media were performed in various systems 0- 
[Toj . In 1989, a non-traditional idea for observing localiza- 
tion was proposed: the transverse localization scheme 0], 
which exploits the equivalence between the Schrodingcr 
equation and the paraxial wave equation for light. In- 
deed, this scheme was used for a clear demonstration of 
Anderson localization in random optical lattices [ll| , for 
the observation of Anderson modes [Hj , and localization 
near an interface [l3|- All of these have dealt with fully 
coherent waves only. In a different domain, the propa- 
gation of partially-incoherent light in random media is a 
subject of considerable interest (e.g., see [lO])- However, 
in the context of localization, the only studies on inco- 
herent light were on enhanced backscattering [2l|, [2^ . 
which is considered a precursor to Anderson localization. 
To the best of our knowledge, strong localization with 
partially-incoherent waves has never been studied. 

Consider the propagation of a partially-spatially inco- 
herent optical beam, linearly polarized, originating from 
a quasimonochromatic continuous-wave source. Such 
a beam can be constructed by sending a laser beam 
through a rotating diffuser (e.g., see [2J]). Its state at 
a given propagation distance z can be described in terms 
of the mutual coherence function uM. 



B{xi,X2,z) = {E* {X2, z,t)E{xi, Z,t))t, 



(1) 



where (• • • }t is the time-average, and E is the stochastic 
field. Instead of B{xi,X2, z), the state of the system can 
be described by an orthonormal set of modes ipj{x,z) 
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(j = 1,2,...), and their modal weights A^, which are 
obtained from the eigenvalue equation [Tsj 



(2) 



i.e., B{xi,X2,z) = J2j >^j'il^jix2,z)-ipjixi,z). 

We analyze linear propagation of such a beam in the 
transverse localization scheme Q, in a waveguide array 



defined by n 



2nQSn{x)^ where ng is a constant 



term, while Sn{x) describes disorder. The propagation of 
the beam along the z axis is governed by the Schrodinger 
equation [ll| 
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1 Sn{x)k 
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(3) 



where k = uquj/c is the wave vector, uj is the temporal 
frequency of the beam carrier, and c is the speed of light. 

In our simulations, we analyze the evolution of Gaus- 
sian input beams: 



B(xi,X2,z = 0) = /oexp[ 



3^1 + X2 
2(7/ 



Xi - X2 



where aj and ac are the spatial and the coherence widths 
of the beam, respectively. In disordered media, the prop- 
agation depends on the particular realization of the ran- 
dom potential Sn(x). Thus, to obtain meaningful results, 
one needs to observe the evolution of the mutual coher- 
ence for many realizations of the disorder, and calculate 
the disorder ensemble average [ll|: X2, z)}^. 

The disordered potential used in our simulations is il- 
lustrated in Fig. [Ija). The index of refraction varies 
randomly (with uniform distribution) between dn — 
and 5ti = 1.2 x 10~^. The width of every rectangular 
potential unit shown in Fig. [Ija) is 2.7 fim, but their 
mutual distances are random: first we fix the leftmost 
rectangle in its position, and then add the adjacent ones 
to the right such that their distance (center to center) 
is between 5 and 9 yum (chosen at random), and so on. 
Such a disordered medium can be created experimen tally 
by using the ultrafast direct laser writing technique |13l |. 

Because the system is linear, the evolution of the beam 
(in a given realization of the disorder) can be described 
in terms of the modes of the system u„ and their propa- 
gation constants which obey 



^ 1 /// N Sn(x)k 

PnU„ = -TrU„{x) H u„ 

Zk no 



(5) 



Every initial coherent wave ipj{x, 0) is projected onto the 



system's modes u. 

ll>j(x, z) = Cj^n^n 



= / dxipj{x,0)ul{x), yielding 



B{xi,X2,z) 



c,.„u„(a;)e*^"^, that is, 

u'^{x2)um{xi)e''^'^" 



E 



(6) 
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FIG. 1: The disordered potential and its eigenmodes. 
(a) Small section of the dimensionless disordered potential 
V{x) — —25n{x){kxoY /no, where k = 27rno/A, A = 514 nm, 
xo = 1 fim, no — 1.45. (b,c) Typical eigenfunctions (Ander- 
son states) of the structure with a smaller (b) and larger (c) 
spatial extent, in a sample of width L = 4 mm. (d) Loga- 
rithmic plot of |«„(x)| from (b) with exponentially decaying 
tails, which are the fingerprint of Anderson localization. 



It is reasonable to conjecture that after sufficiently long 
propagation z, the beam will dephase and attain the fol- 
lowing form of the mutual coherence: 



B(xi,X2, z 



^Aj|cj,„p<(a;2)w„(a;i). (7) 



From Eqs. ^ and ([7]) we clearly see that, if the medium 
is infinitely broad, and if the potential is fully random, 
then, since all the eigenmodes u„ are Anderson localized, 
any initial finite-width partially-incoherent beam will not 
diffuse during propagation despite its incoherence. One 
can expect that, after sufficiently long propagation, the 
tails of the incoherent beam will become exponentially 
decaying: B{x,x,z — >■ oo) oc exp(— 7|a;|), with 7 cor- 
responding to the excited eigenmode u„ which has the 
slowest decay. However, in reality, all samples have fi- 
nite transverse (L) and longitudinal {Z) size; hence, we 
explore the finite-size effects. 

First we calculate the eigenvectors/eigenvalues of the 
potential 5n{x). We use the LAPACK implementation 
of the MRRR algorithm [23|, with Dirichlet boundary 
conditions m„ = at the boundaries of the medium 
{x = ±L/2). Typical eigenfunctions are shown in Figs. 
[IJb)-(d). Up to some critical value of n, the spatial extent 
of the eigenfunctions u„ is smaller than L [see Fig. [IJb)]. 
whereas above this value they extend to one or both of 
the boundaries [sec Fig. [Ijc)]. The localized eigenfunc- 
tions are a consequence of Anderson localization, which 
can be seen from Fig. [T](d), where we depict |m„(x)| on 
a logarithmic scale. Clearly, the amplitude of Un{x) de- 
cays exponentially, log |u„(a;)| oc exp(7„a:), which is a fin- 
gerprint of Anderson localization in this (l-l-l)D system. 
From these eigenmodes, we calculate their Lyapunov ex- 
ponents as follows: We form a set {|un(a^max)|} of local 
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FIG. 2: (a) Lyapunov exponents 7„ extracted from the eigen- 
functions Un (solid blue line) and via the transfer matrix ap- 
proach (black squares). Shown are averages over 40 different 
realizations of the random potential; the standard deviation 
is smaller than the thickness of the line, (b) The Thouless 
(solid blue line) and the Heisenberg (dashed red line) times 
(propagation distances) vs. n for the system with L — 4 mm. 
Horizontal line depicts Z — 10 cm. See text for details. 



maximaof |it„(a;)|, and then fit — 7„|a;— xd+fe to the set of 
points log |Mn(a;inax)|, to obtain 7„ for every m„. For the 
lowest eigenstates, which have too few local extrema, |u„| 
was fitted. The Lyapunov exponents, averaged over 40 
different realizations of the random potential, are shown 
in Fig. HJa) as a function of n; the standard deviation 
is as small as the thickness of the line implying that the 
Lyapunov exponents have the seff-averaging property. In 
order to underpin our results, we have also calculated the 
Lyapunov exponents via the transfer matrix approach 
[l9|, i.e., for an essentially infinite sample, and obtained 
nondistinguishable results [black squares in Fig. [DJa)]. 
Evidently, the overall trend is such that the Lyapunov 
exponents decrease with increasing n (decreasing /3„), 
except for a small hump, which is a consequence of rel- 
atively small fluctuations in distances between adjacent 
potential peaks [Fig. [IJa)] (if we allow larger fluctuations 
the hump disappears). 

In a realistic experimental situation with finite-size 
samples, a partially-incoherent initial wavepacket is likely 
to excite both types of modes: modes that are Anderson 
localized on a scale smaller than L, and modes extending 
to the boundary. In order to investigate the finite-size ef- 
fects, we first evolve (numerically) an initial wavepacket 
with absorbing boundary conditions [Figure |3] (a)], and 
then repeat the simulation with reflecting boundary con- 
ditions [Figure |3] (b)] (both can be realized in experi- 
ments). The initial size of the beam is aj ~ 10 /im, and 
the degree of coherence varies: ac = 1,3, and 5 /im; 
L = 4 mm. In our simulations, absorbing boundary con- 
ditions are included as an imaginary index of refraction, 
QSriahs > 0, which is present only close to the sample 
edges: Sn^hs = for |.t| < 0.96f . 

From the simulations with absorbing edges we find 
that, after sufficiently long propagation, the intensity 
structure has exponentially decaying tails: {I)d oc 
exp(— 71x1). All graphs have approximately the same 
value for the slope, 7 = 3.2 x 10"'^ fim~^. This is in 
accordance with our analysis of eigenstates. Namely, the 
part of the beam exciting eigenstates which touch (or are 
in the very vicinity of) the edges, gets absorbed during 
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FIG. 3: Intensities of the incoherent beam after Z — 10 cm 
of propagation, for different system parameters, averaged over 
40 realizations of the disorder. (a,b) Beam intensity for ab- 
sorbing (a) and reflecting (b) boundary conditions for beams 
with a I = 10 /im, and ac = 1 /im (top black line), 3 /tm (mid- 
dle blue line), and 5 /tm (bottom red line); L = 4 mm. (c) 
Beam intensity for absorbing boundary conditions and three 
values of L: 2, 4, and 8 mm (top to bottom) [aj — 10 /im, 
ac = 1 /im]. (d) The intensity for one specific instantaneous 
initial realization of the incoherent field (lower blue line), and 
the time-averaged intensity (upper black line) [aj = 10 /im, 
ac ~ 1 /tm, L = 8 mm, absorbing boundary conditions]. 



propagation. The remaining part of the beam excites 
only exponentially decaying eigenstates, which yields the 
intensity plotted in Fig. ^a.). The value of 7 found 
from graphs in Fig. El^a) correspond to the slowest de- 
cay rate 7„ of the eigenstate it„ which is excited, and 
which does not overlap with the absorbing boundary, i.e., 
{un\Snahs{x)\un) ~ (ouc should take into account factor 
2 because /= {\E\^)t)- 

At this point it is worthy to consider the Thouless (Zt) 
and Heisenberg (Zh) times (propagation distances) in 
our system. The former tells us the average time it takes 
for a particle (with some energy) to diffuse across the 
sample, and the latter is the longest time the particle 
can travel inside the sample without visiting the same 
region twice (e.g., see [ij] and Refs. therein). In our 
case Zt corresponds to the inverse linewidth (l/A/3„) of 
the eigenstates with absorbing boundary conditions. In 
calculating Zxin) we have averaged l/A/3„ over 20 ad- 
jacent eigenstates, and then over 40 realizations of the 
potential. In the same fashion we have calculated the 
average inverse level spacing which yields the Heisen- 
berg time Znin). Figure [DJb) shows Zt and Zh vs. 
n. Evidently the Thouless time is effectively infinite for 
n < 1200 indicating Anderson localization in our finite 
sample, whereas for greater n values the two times are 
on comparable scales. From Fig. Efa), we conclude that 
the more coherent the light is - the stronger the localiza- 
tion effect is, in a sense that less power is absorbed by 
the walls. For the samples used in our simulations the 
rate of decaying tails are identical for all values of ac, 
because in every case the highest of the localized modes 
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FIG. 4: Spatial power spectrum of the incoherent beam after 
Z = 10 cm of propagation in a disordered medium with ab- 
sorbing (thin dashed blue line) and reflecting (solid red line) 
edges, for an initial beam (thick solid black line) defined by 
ai = 10 nm, Gc = 1 ^m. All plots are averages over 40 
realizations of the disorder. 



was excited. 

It is interesting to note that the values 7 obtained after 
transport simulations (expansion) depend on the size of 
the sample L, and scale approximately as 7 oc L~^. We 
checked this trend for L = 2, 4, and 8 mm. see Fig. ^c). 
The observed scaling is found to be in agreement with 
the studies conducted in 



25|, where it was shown that 



the logarithm of the dimensionless resistance scales as 
the length in case of localization. 

Figure [4] shows a typical change in the spatial power 
spectrum (closely related to the spatial coherence) for 
propagation dynamics with different boundary condi- 
tions. In the case of absorbing edges, there is a clear 
cut off in fc-spacc above which the modes are no longer 
localized on a scale smaller than L. The modes u„ com- 
prised of plane waves (via Fourier transform) with k val- 
ues above the cut-off evidently spread to the absorbing 
boundaries of the medium. This process also increases 
the overall coherence of the beam, because it has fewer 
and fewer modes as it propagates. In contrast, for re- 
flecting boundary conditions, the final power spectrum 
has the same overall shape as the initial spectrum, with 
fluctuations on top of the average shape. 

A partially-spatially incoherent (quasimonochromatic) 
beam can be thought of as a train of fully coherent fields 
with intricate field structures, which are replaced on a 
time scale ~ Aic (the temporal coherence time); quasi- 
monochromatic means that Aic 3> Z/c. In experiments, 
one instantaneous realization of the incoherent field can 
be observed by simply stopping the rotation of the dif- 



fuser [231 . Interestingly, we find that one instantaneous 
field realization will asymptotically have exponentially- 
decaying tails similar to those of the time-averaged inten- 
sity, as illustrated in Figure (both graphs are aver- 
aged over disorder). Thus, the localization of a partially- 
incoherent field can be thought of as the localization of 
many instantaneous coherent speckled structures, whose 
ensemble average produces localized mutual coherence. 
The exponential decay of the asymptotic intensity evi- 
dently depends both on the degree of coherence and the 
spatial structure of the input field. 

We can extend our conclusions to hold for spatially and 
temporally incoherent beams of light (Ate < Z /c)^ such 
as the one constructed from an incandescent light bulb 
in Ref. [2^. In this case, a space-time mutual coher- 
ence function describing the beam can be transformed 
into the space-frequency domain [ll, [23|. Because our 
results hold for every frequency component separately 
(with somewhat different numerical values because dif- 
ferent frequencies experience different diffraction), due 
to linear superposition, it is evident that they also hold 
for the whole spatially and temporally incoherent beam. 

Before closing, we note that our results correspond 
to partially condensed noninteracting Bose-Einstcin con- 
densates (BECs) in ID disordered potentials. As an in- 
teresting consequence of dimensionality, we note that our 
results also pertain to the strongly interacting ID Bose 
gases [2^, and noninteracting ID Fermi gases, which are 
modeled by identical equations [2^ . To mimic partially- 
condensed interacting BECs [13] , one should include non- 
linearity, which we leave for further studies. 

In conclusion, we have predicted the properties of 
Anderson localization of partially spatially-incoherent 
beams in disordered linear (l-l-l)D photonic structures. 
We conclude that more incoherent light will diffusively 
spread more through the random medium than coher- 
ent light, however, the incoherent wavepacket will dis- 
play exponentially decaying tails after sufficiently long 
propagation. We have discussed the finite size effects, 
and extended our conclusions for spatially and tempo- 
rally incoherent light beams. 
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